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Рассмотрена осесимметричная контактная задача теории упругости о кручении круговым штампом трансвер- 
сально-изотропного полупространства с неоднородным покрытием. С помощью техники интег! ральных пре- 
образований решение задачи сведено к решению интегрального уравнения. Построено явное выражение 
трансформанты ядра интегрального уравнения для однородного трансверсально-изотропного покрытия. 
Предложена схема численного построения трансформанты ядра в случае произвольного закона изменения 
модуля сдвига по Г лубине покрытия. Используя аппроксимацию трансформанты ядра некоторыми специаль- 
НЫМИ функциями, построено двухсторонне асимптотически точное приближённое аналитическое решение 
задачи. Проанализирован частный случай изменения модулей сдвига в трансверсально-изотропном покры- 
ТИИ, ПОЗВОЛЯЮЩИЙ построить явное замкнутое решение задачи. Численно исследуется модельная задача о 
кручении биоматериала. Считается, что покрытие представляет собой хрящевую ткань, а свойства подложки 
соответствуют кортикальной КОСТИ. 
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Введение. Первые работы по моделированию задачи о кручении упругого полупространства 
круглым штампом относятся к середине прошлого века. Задача о кручении однородного упругого 
полупространства круглым штампом впервые была сформулирована и решена в динамической 
постановке Рейснером и Сагочи [1]. В работе Снеддона [2] данная задача, с использованием тех- 
ники интегральных преобразований, сведена к решению интегрального уравнения. 

Грилицкий В. Д. [3] построил решение задачи о кручении круглым штампом изотропной 
двухслойной среды и ортотропного упругого слоя в виде ряда по степеням б`!, где б — отноше- 
ние толщины первого слоя к радиусу штампа. Каз$! М. К. [4] построил решение задачи о круче- 
нии неоднородного изотропного полупространства в случае степенного изменения модуля сдвига 
по глубине. Решение задачи для изотропного полупространства с неоднородным покрытием при 
произвольным законе изменения модуля сдвига по глубине получено Айзиковичем С. М. [5]. 
Егдцуеп М. Е. рассмотрел задачу о кручении трансверсально-изотропного упругого полупростран- 
ства с однородным покрытием [6]. 

Задача активно изучается по настоящее время, в частности, \№Мапд \.-5. использовал ку- 
сочно-линейную аппроксимацию трансформанты ядра для сведения задачи к решению инте- 
грального уравнения [7]. В работе [8] решена задача о кручении трансверсально-изотропного 
пьезоэлектрического однородного полупространства. 


1 Работа выполнена при поддержке Минобрнауки РФ (ГК № 11.519.11.3028, соглашения №14.В37.21.1632, 14.В37.21.1131, 
14.132.21.1693) и РФФИ (12-07-00639_а). 
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В настоящей работе построено решение контактной задачи о кручении трансверсально- 
изотропного упругого полупространства с неоднородным трансверсально-изотропным покрытием 
при произвольном изменении упругих свойств по глубине покрытия. 

Изменение механических свойств по одной из координат и анизотропия упругих свойств 
характерны для многих тел, как в природе, так и в технике. В технике неоднородность и анизо- 
тропия связаны с условиями создания и эксплуатации деталей машин. Большинство биологиче- 
ских объектов являются неоднородными и анизотропными, что свидетельствует об актуальности 
построения математических моделей, учитывающих одновременно и неоднородность, и анизо- 
тропию. 

Постановка задачи. Недеформируемый круглый штамп с плоским основанием жёстко сцеплен с 
верхней гранью Г упругого неоднородного трансверсально-изотропного полупространства ©. С 
полупространством связана цилиндрическая система координат г,ф,2. Штамп контактирует с 


полупространством по поверхности 2 =0, г<а. 





Рис. 1. Задача о кручении 


К штампу приложен крутящий момент М, ось которого совпадает с осью 2 (рис. 1). Под 
действием этого момента штамп повернётся относительно оси 2 на угол =, вызвав деформацию 
кручения ©. Модули сдвига С„„ Су. полупространства с глубиной изменяются по законам: 


Е (2) -Н<2<0 Е (2) -Н<2<0 


гф фр 


< = 1 
“ В.Е. (-Н) —<2<-Н, “ [В,-Ё,(-Н) —<2<-Н, (1) 


где Го, »› — некоторые функции, определяющие закон изменения модулей сдвига в покрытии, 
причём Г» — непрерывная, а №, — непрерывно дифференцируемая; параметры Ви, В> задают 
скачок упругих свойств на границе покрытие-подложка. 
Вне штампа поверхность Г не нагружена. При сделанных предположениях граничные 
условия задачи имеют вид: 


Т, =0, Г>а 


2=0: в, =Т, =0, (2) 


Учитывая, что слой подвержен лишь деформации скручивания, ненулевым является лишь смеще- 
ние вдоль оси ф: 


и, =гЕ, гза. 


и. =И, =0, и, =и(г,2). 
Напряжения, в соответствии с законом Гука, примут вид: 


ди ди, и 
9, =0,=9, =Т, =О,Т., =б., (2) и То = бу т". (3) 
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Введём обозначения: 


т —Н<27<0 и —Н<2<0 
т, = и, = : 
} т 2<-Н, ы в 2<-Н 


При г > о и 2 > —ю напряжения исчезают. Считаем, что перемещения и напряжения сопряга- 
ются на границе изменения закона неоднородности: 
. 1 2 1 2 
2=-Н: тет и® =и@. (4) 


Ф2 Ф Ф 

Требуется определить закон распределения контактных касательных напряжений под штампом 
т ет. (г), г<а. 

Сведение задачи к решению интегрального уравнения. Уравнения равновесия при данной 

постановке принимают вид: 














р би 
+ ИО |+ — + +С =0, -Н<2<0, 5 
.| дг дот 4 02 ть (5) 
10и® 0и® 1 218) 
©. (-Н и и@ +6 (НН) =0, 2<-Н. 6 
и |. дг дг? [2 Ф Че ) 922 ( ) 


Для сведения задачи к решению интегрального уравнения используем интегральное преобразо- 
вание Ханкеля: 








4 (г.г) = [0 (у,2)(су)уа, 1=Ъ2. 7 
Обозначим 
С. (2 С, (2 
еб, 5-Е) о 


Подставив (7), (8) в (5) и (6), и приравнивая к нулю в полученных соотношениях подынтеграль- 
ные выражения, получим систему обыкновенных дифференциальных уравнений относительно 


функций И, (\,2), 1=1,2: 
(у, 2) +у(2) Ц, (у,2)-У?к(2 Ц, (у,2)=0, -Н<2<50, (9) 
И (у,2)- ук (-Н)Ц, (у,2)=0, -®<2<-Н. (10) 
Общее решение уравнения (10) имеет вид: 
и, (\,2) = В, (у) С® + в, (уе СВ. 
Из условия ограниченности перемещений при 2 -> —ю имеем В, (у) =0 , тогда 
О, (\,2) = В, (у)е*\С#) : (11) 
Условия сопряжения (4) можно записать в виде: 
и, (у,-Н) = И, (у,-Н), Ц, (у,-Н ) = ВМ; (у,-Н). 
Или, используя явный вид функции И, (у,2 (11), получим: 


(%,2) 
те - УКС (12) 


И, (у, 
Из (3) при 2 = 0и (7) следует: 
И; (\,0) =Т (\)/6. (0). (13) 
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Здесь Гу) связана с т, (г) преобразованием Ханкеля: 


)= [т, (2), (ру)рар. (14) 

0 

Введём вспомогательную функцию (/" — ): 
Ш: (у,2) = УЦ, (у,2)Т^ (У) 6» (0)6. (0). (15) 
Ввиду однородности, уравнение (9) относительно функции И’ (у,2) не изменится. Условия (12), 


(13) относительно И’ примут вид: 
и” у, -Н) 
в ВЫ = Ву. к (- (16) 
о 
0" (у,0) =Ук( (17) 


Далее Ц; (у,0) будем обозначать / У. — используя условие (2), можно записать инте- 


Гральное уравнение нашей задачи: 


© 


о И, (у, 0) 7, (гу) уау = 
м-р 


Используя (14), уравнение можно представить в виде: 
| (2)в] Е (У)-7 (гу) 1 (ру) дуар =г= 6, (0)б., (0), г<а. 
0 0 


Сделаем замену переменных: 


Ну=и, ^=Н/а, г'’=г/а, р'=р/а,Ё(и)=Ё(и/Н),т(р’)=т, (р’а). 
Далее штрихи опускаем, тогда интегральное уравнение задачи примет вид: 


но) (оур[и (и) , (шел), (ир^”" )аиар = №= 6, (0). (0), г<1. (18) 


Таким образом, для построения интегрального уравнения задачи необходимо построить функцию 


Ч 








Не 1, (гу)а4у =г=, г<а. 
)° 


И (у,0) из краевой двухточечной задачи (9), (16), (17). 


Построение трансформанты ядра интегрального уравнения. Для некоторых частных слу- 
чаев изменения модуля сдвига в покрытии удаётся получить явное выражение для трансформан- 
ты ядра [(и). Например, для покрытия, в котором модули сдвига постоянны: 


66, -Н<2<0 6, -Н<2<0 
бу (2)= т ( )= Е 


гф! фг! 
В этом случае решение уравнения (9) можно построить в явном виде, и трансформанта ядра при- 


нимает вид: 
122% 665-655, 6с 
К(и)= ин | Н= с к # б=\|-с`›б. = 
1-не Сб, +60, С 


Заметим, что для многослойной среды уравнение (9) распадается на № независимых дифференци- 
альных уравнений с постоянными коэффициентами, для каждого из которых известно явное ре- 
шение. Для построения значения трансформанты ядра в некоторой точке достаточно решить си- 
стему линейных алгебраических уравнений. 

Для произвольного закона изменения модуля сдвига по глубине значение трансформанты 
ядра может быть получено численно. Введём обозначения: 


0:= М (2); =): 
28 


-®<2<-Н, —ю<2<-Н. 








Вестник ДГТУ. 2013. № 1—2 (70—71) 








Уравнение (9) перепишем в виде линейной системы обыкновенных дифференциальных уравнений 
первого порядка: 


0; =0, 
' 2 5 Н 2 < 0 Г 
в» =-у(2)6, +У'к(2)6:, 
или в матричной форме: 
О, 0 1 

46 А, -Н<2<0, А-|, И 
[2 уж(2) —*(2) и 
Краевые условия (16), (17) принимают вид: 


в," |„_и =В.\\К(-Н), @, |. =У\К(0). (19) 
Решение краевой двухточечной задачи при фиксированном у построим, используя метод модели- 
рующих функций [9], предложенный для построения трансформант ядер динамических контакт- 
ных задач. Ищем &(у,2) в виде: 
@(у,2) = В, (у)а (у,2)е\. 
Вектор а (\у,2) определяется из задачи Коши при фиксированном у: 
а 
[674 
а, |р > В, ‚ а | — Ук (-Н) . 


Затем, используя второе условие (19), находим В, (У) итем самым ‹, (у,0): 
В, (у) = уа — к(0), 


Кин) =Е (у) = (у,0) = в (у,0) = ук (@)а, (+ б)аз (у,0). 


Аналогично работе [10] можно показать, что о ядра и р: свойствами: 


=Аа-уа, -Н<2<0, 





ЕЕ: т (и) =1. (20) 


Построение решения задачи. Будем говорить, что функция Хи) принадлежит классу Пл, где 
№ — натуральное число, если она представима в виде: 
м 1+ А 
(и) = —— 
( ) | и? Е В? / 
Для трансформанты ядра [, (и) ЕП» можно построить [11] решение интегрального уравнения 
(18) в форме: 


А,В ЕС, А*А, В, +В, при! =]. 


4 1 
т(/) =—=\/б (0)С.„ (0 С, (г. АЛ), (21) 
рее с2Ар 

где постоянные С, определяются из системы линейных алгебраических уравнений: 

м А В 1+В Л" 

С Ра К 3 К =0, К=1,2,...,М. 

>: 29 э г, (0) В?^? 
Здесь использованы обозначения: 
$И Аг ГНА Аг | св АЕ 


Китай) [+ -п) 
А’). 


р(А,В) = (АА+ВЗНА)/(В? 





2(г,А)= 
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Таким образом, для построения решения исходной контактной задачи достаточно аппроксимиро- 
вать трансформанту ядра интегрального уравнения функцией класса Пл. Для построения аппрок- 
симации можно использовать алгоритм, описанный в работе [11], где показано, что решение (21) 
является асимптотически точным при Л-0 и Л-›09. 

Решение задачи для частного случая изменения модулей сдвига. Пусть для модулей 
сдвига выполнены равенства: 








бь(2)б., (2) =с* = сопз, В, =В, =1. (22) 
Тогда первое равенство (8) можно переписать в виде: 
к(2) =б„(2)/с. (23) 
Учитывая (23), уравнение (9) принимает вид: 
6. (2), ИР(у,2) у? 
- 6? (Е 0! (у,2) + 6. Е еб (2), (у,2)=0, -Н<250. (24) 
Общее решение уравнения (24) имеет вид: 
уз) = ес Е", (25) 


где Е (2) = [6 (2). 


Продифференцируем ( (у,2) по 2: 


а . 2)\/С -Е(2)у/с 
ОЕ (м2) = 6, (2) [6 (уе”®" 6, (де |. 


[674 

Тогда из краевого условия (16), учитывая (23), следует: 

Е. (у) =0. (26) 
Из краевого условия (17) и равенств (23), (26) получаем следующее: 

с, (\)е"® =1. (27) 
Из (25), (26), (27) следует, что для любого у выполнено тождество 
и: (у,0) =1. 

Тогда 

ГВ (и) =1. 


Таким образом, если для модулей сдвига трансверсально-изотропного полупространства с неод- 

нородным трансверсально-изотропным покрытием выполнены равенства (22), то трансформанта 

ядра интегрального уравнения тождественно равна единице. Следовательно, решение задачи 

совпадает с решением для однородного полупространства с модулем сдвига равным си имеет вид: 
Г 

Численные примеры. Рассмотрим трансверсально-изотропную упругую подложку с модулями 

сдвигаб,(-Н)=5, б..(-Н)=6 и пусть модули сдвига на поверхности покрытия равны 





т(г) = пес 


Сб (0)=С., (0) =0,608 . Такие значения модулей сдвига соответствуют свойствам хряща (покры- 


тие) и кортикальной кости (подложка). Предположим, модули сдвига покрытия изменяются по 
одному из следующих законов: 
1. (2) =5,03- 4,422е°4" , ,(2)=6,037 -5,428е°*"" ; 
2. 1.(2)=0,608 - 4,3922/Н , #,,(2)=0,608 —5,3922/Н ; 
3. #.(2)=#, (2) = 0,608; 
4. Е. (2)=Е„ (2)=0,118+0,489е°4* . 
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Эти законы описывают трансверсально-изотропные покрытия, модули сдвига которых возрастают 
экспоненциально (закон 1) и линейно (закон 2); однородное изотропное покрытие (закон 3) и 
изотропное покрытие (закон 4), упругие свойства которого экспоненциально убывают, а на гра- 
нице покрытие-подложка наблюдается скачок модулей сдвига Си, Сь, в 41 и 49 раз соответствен- 


но (В:=41, В›=49). Графики изменения модулей сдвига по глубине для законов 1—4 изображены 
на рисунке 2. 


0 0 
я 28. 
Н 1 Н 1 
-0.5 = 0.5 
-1 -1 


0.1 0.6 Ср: э 6 0.1 0.6 Ср а) 6 


Рис. 2. Изменение модулей сдвига по глубине в материалах 1—4 


Трансформанты ядер, соответствующие покрытиям 1—4, изображены на рисунке 3. По- 
строены аппроксимации трансформант ядер функциями класса Пл, их погрешности равны 0,8 %, 
0,82 %, 0,64 % и 0,49 % соответственно. Под погрешностью аппроксимации понимаем следующее 
выражение: 


А, = тах [1 (и)/ (и)-1|х100%. 


-2 -1 0 
10 10 10 10 10 10 и 10 
Рис. 3. Трансформанты ядер для материалов 1—4 


Введём величину: 
ЕР а 
где ^=Н/а — геометрический параметр задачи, характеризующий относительную толщину по- 
крытия; т(^, /) — контактные напряжения под штампом для рассматриваемой неоднородной 


среды; Тот (г) — контактные напряжения однородной изотропной среды с модулем сдвига, 


равным 6 (-Н)6.. (-Н). 





Величина Т,. характеризует распределение контактных касательных напряжений под 
штампом для неоднородной среды по сравнению с однородной средой. 
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Графики величины Т.. для материалов 1—4 в случае относительно толстого покрытия 


ге] 
(Л=100), тонкого покрытия (Л=1/100) и покрытия, сравнимого по размерам с размерами штампа 
(Л=1), изображены на рисунке 4. Из графиков видно, что неоднородность покрытия оказывает 
существенное влияние на решение контактной задачи, особенно для покрытий малой толщины. 
Наиболее существенное перераспределение контактных напряжений по сравнению со случаем 
однородного полупространства наблюдается на краю штампа (1). 
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Рис. 4. Контактные напряжения для законов 1—4 


Заключение. Построено приближённое аналитическое решение контактной задачи о кручении 
трансверсально-изотропного упругого полупространства с неоднородным трансверсально- 
изотропным покрытием, эффективное для любых относительных толщин покрытий (значений па- 
раметра Л). Показано, что в случае, когда модули сдвига трансверсально-изотропного полупро- 
странства изменяются в соответствии с условием (22), распределение контактных напряжений 
под штампом совпадает со случаем однородного изотропного полупространства. 

На основании высокой точности аппроксимации трансформант ядер интегральных уравне- 
ний (менее 0,8 %) можно утверждать о малой погрешности построенных решений контактной за- 
дачи даже для покрытий сложной структуры (материал 4). 

Построенное в работе решение позволяет исследовать одновременное влияние анизотро- 
пии и неоднородности покрытия на механические характеристики контактного взаимодействия. 
Результаты актуальны для широкого класса современных полимерных, биологических и компози- 
ционных материалов. 
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